Optimal stochastic power flow using enhanced multi-objective mayfly algorithm

For the classical multi-objective optimal power flow (MOOPF) problem, only traditional thermal power generators are used in power systems. However, there is an increasing interest in renewable energy sources and the MOOPF problem using wind and solar energy has been raised to replace part of the thermal generators in the system with wind turbines and solar photovoltaics (PV) generators. The optimization objectives of MOOPF with renewable energy sources vary with the study case. They are mainly a combination of 2–4 objectives from fuel cost, emissions, power loss and voltage deviation (VD). In addition, reasonable prediction of renewable power is a major difficulty due to the discontinuous, disordered and unstable nature of renewable energy. In this paper, the Weibull probability distribution function (PDF) and lognormal PDF are applied to evaluate the available wind and available solar power, respectively. In this paper, an enhanced multi-objective mayfly algorithm (NSMA-SF) based on non-dominated sorting and the superiority of feasible solutions is implemented to tackle the MOOPF problem with wind and solar energy. The algorithm NSMA-SF is applied to the modified IEEE-30 and standard IEEE-57 bus test systems. The simulation results are analyzed and compared with the recently reported MOOPF results.


Introduction
The optimal power flow (OPF) problem was first solved by an optimization method with exact second order derivatives in 1984 [1], and its definition is network constrained economic OPF was provided in Ref. [2].The purpose of the classical OPF problem is to optimize a certain objective or multi-objective of minimizing fuel cost, emission, real power loss, VD, etc.While minimizing single or multi-objectives, feasible control variables are calculated to satisfy various equality and inequality constraints.However, with the large-scale use of generators in thermal power plants, large amounts of harmful emissions are released, thus causing great harm to the environment [3].Emissions are also universally considered as a necessary goal to reduce environmental pollution.Hence, MOOPF problem consists of fuel cost, emission, power loss and VD in this study.Due to the highly nonlinear nature of the MOOPF and the • An improved multi-objective mayfly algorithm (NSMA-SF) is proposed.In NSMA-SF, three strategies are incorporated to enhance the capabilities of the basic mayfly algorithm.• The enhanced mayfly algorithm to solve the MOOPF problem with wind and solar energy on two different bus systems, and the robustness, superiority and efficiency of the algorithm are demonstrated.• The experiment results show that the proposed NSMA-SF can be an effective for MOOPF problem.
The rest of study is organized as follows.The mathematical model of MOOPF with stochastic renewable power is introduced in Section 2. Mathematical models of wind and solar power generators are described in Section 3. In Section 4, the multi-objective algorithm NSMA-SF is proposed for MOOPF with wind and solar energy.Section 5 shows the simulation results and comparative studies.The conclusion and future work are presented in Section 6.

The MOOPF problem with renewable energy sources model
Mathematically, the MOOPF consists mainly of the selected optimization objective f i (x,y), control and state variables (x and y) and constraints (g(x, y) and h(x, y)), which can be formulated as Eq.(1).

Control variables and state variables
The MOOPF problem with wind and solar energy includes the control (or independent) variables, which is composed of the active power, voltage magnitude, and shunt compensators.The control variable x is represented in the vector form as Eq. ( 2): x = [ P GB2 , ..., P GBNG , V GB1 , ..., V GBNG , Q C1 , ..., where P GBi is the active power of the generator at the i-th bus (except the slack generator P GB1 ), V GBi denotes the voltage magnitude, Q Ck represents the shunt compensation at k-th bus.The state variable y represents a quantitative description of the actual state of the entire power system network.The state (or dependent) variable y is defined in vector form as Eq. ( 3): where V LBi denotes the actual voltage at i-th bus.The reactive power of the generator is written in Q G , S l represents the loading of the transmission line.

Constraints
Constraints are classified as equality constraints and inequality constraints.Specifically, equality constraints serve primarily to regulate the equilibrium of power input and output, whereas operational constraints give rise to the inequality constraints.

Equality constraints
In a power network, the actual active power value of the generator is as large as the reactive power output.The equality constraints formulated as Eqs.( 4) and (5): where θ Bij denotes the angle difference between the voltage phasor values on the i-th and j-th bus, NB indicates the number of total buses, P DBi is the power demanded of the generator, Q DBi represents the reactive power demanded of the generator, G Bij and B Bij are identified as the transfer conductance and susceptance values between i-th and j-th bus.

Inequality constraints
Inequality constraints are based on inequality constraints caused by equipment operating constraints can be defined as Eqs.( 6)-( 12).

Generation cost
In the MOOPF, the minimization of the generation cost is the main objective.The generation cost is defined as the following formulas ( 13) and ( 14): where E Gi denotes the energy cost of valve point loading effect in traditional thermal power generators, a i , b i and c i denote the corresponding generator cost coefficients.d i and e i indicate the valve load effect coefficients.The specific values of each coefficients can be presented in Table 1.

Emission
With the increasingly stringent environmental protection requirements, reducing emissions is a reasonable optimization objective.Conventional thermal power generation will emit a large amount of greenhouse and harmful gases.The emission is calculated as Eq. ( 15), the specific values of emission coefficients are presented in Table 2.

Real power loss
In the transmission of a real power system, the power loss (or outages) in the system cannot be avoided, and power loss can be calculated as Eq. ( 16):

Voltage deviation (VD)
For a power system, the VD is a significant indicator to detect the voltage quality.The total VD is formulated as Eq. ( 17):

Case studies
In this paper, 10 different multi-objective case studies are conducted in the modified IEEE-30 and the classic IEEE-57 bus system.The multi-objective of the MOOPF problem is composed of two to four objectives.Table 3 summarizes the study cases.
3. Stochastic wind and solar power model.Section 3 is focused to an efficient approach for modeling the renewable energy generators.Replacing the thermal generators with wind turbines and solar PV farms and calculating the cost of two generators.Fig. 1 is a wiring diagram of the modified IEEE-30 bus test system.

Wind power model
In a multitude of simulation experiments, the probability of general wind speed is divided into a skewed normal distribution.With in-depth research, the two-parameter Weibull distribution has been as the better function suitable for the actual wind speed distribution.The Weibull PDF is widely used in wind energy resource assessment and economic evaluation of a wind farm.

Stochastic wind speed probability distribution function
Referring the two-parameter Weibull distribution, the wind speed PDF is given as Eq. ( 18): where v denotes the actual wind speed, λ and k denote the scale and shape parameter.Wind speed distribution for wind farms 1 and 2 are shown in Figs. 2 and 3.
In the 30-bus system, two wind generators with different output powers are used.The actual output power (P w (v)) of each turbine is calculated as Eq. ( 19)： (19) where v in denotes the cut-in wind speed of the turbine, v r expresses the rated wind speed, and v out indicates the cut-out wind speed.P wr indicates the rated power output of the turbine.The specific values of the correlation coefficients are in Table 4.
The discontinuity, disorder and instability of the wind resource results in variable electrical energy from wind turbines.According to the different cases of wind speed, the probability calculation of three different states is given by Eqs. ( 20)- (22):

Direct cost evaluation of wind power
Unlike thermal generators, there are no fuel costs for wind generators, but there are costs to purchase wind turbines, which are called direct costs.The direct cost of a wind generator (C Dw ) is computed as Eq. ( 23):   J. Zhu et al.
where g i expresses the direct cost coefficient, P ws,i denotes the scheduled power and Nw indicates the number of wind generators (Nw = 2).The values of the related parameters are provided in Table 4.

Wind power cost evaluation with uncertainties
Due to the discontinuous, disordered and unstable nature of wind energy, wind turbines do not generate electricity continuously, but intermittently, and operators must also consider some other costs.The additional costs that arise in two cases when wind turbines generate excess and deficit power, namely the reserve cost (C Rw ) and the penalty cost (C Pw ).The calculation formulas for C Rw and C Pw of a wind power plant are displayed in the following formulas (24) and (25): where K Rw,i , C Rw,i are the reserve and penalty cost coefficients, respectively.

Solar PV power model
Due to the great influence of factors such as the intensity of solar radiation and temperature, the power generation time of PV generators has certain fluctuations and uncertainties.Therefore, in the power grid, energy storage systems will be used for complementary power generation to alleviate the uncertainty of the output of PV generators.

Calculation of solar PV power probabilities
Some studies have pointed out that in a short period (within one day), the illumination of the sun can be approximately regarded as a Beta distribution [19].Under the premise that the influence of seasonal changes on light intensity is not considered, a large number of papers show that the log-normal PDF can be applied to fit the probability distribution of solar radiation intensity.The log-normal PDF is formulated as Eq. ( 26):

Table 4
Parameters of wind power generator plants.
where G, σ and μ represent the solar irradiance, standard deviation and mean value, respectively.The values of the coefficients are provided in Table 5.The solar irradiance distribution for solar PV is shown in Fig. 4. In the PV generators, according to Eq. ( 27), solar irradiance can affect energy conversion with the equation: where P sr denotes the rated output power of the solar PV generator.Fig. 5 is the real power distribution of a solar PV generator.

Direct cost evaluation of solar PV power
According to Ref. [20], the direct cost formula for a solar generator (C Ds ) is displayed in the following formula (28): where P ss,j is the scheduled power.h j represents the direct cost coefficient.Ns denote the number of solar PV power plants, in this study Ns = 1.

Solar PV power cost evaluation with uncertainties
Similar to wind power generators, storage costs and penalty costs need to be calculated in the solar PV generator.According to Ref. [20], the service cost (C Rs ) and penalty cost (C Ps ) of solar PV generators are calculated as Eqs.( 29) and ( 30): where K Rs , K Rs denote the reserve and penalty cost coefficient.Specific coefficient values are provided in Table 5.

Modified multi-objective Mayfly algorithm
This section will discuss how to solve the MOOPF problem with NSMA-SF.First, the classical MA algorithm is briefly introduced, and an improved MA algorithm (NSMA-SF) is proposed based on it.In NSMA-SF, nonlinear decreasing inertia weights are introduced, and the fast non-dominated sorting in the NSGA-II is referred to.Furthermore, an efficient constraint processing technique, called the superiority of the feasible solutions (SF), is introduced.

Standard mayfly algorithm (MA)
The MA is a new optimization algorithm with an algorithmic mechanism inspired by mayfly populations [21].Different update behaviors are performed by male and female mayflies.The real world has numerous multi-objective problems, and each objective is often directly conflicting, which cannot be solved by MA algorithms.When the MA algorithm was used to solve the multi-objective optimization problem, the results show that the optimization effect is ineffective.In this case, a better solution can be obtained by compromising each sub-objective.

Multi-objective mayfly algorithm
In the NSMA, as with the standard MA, the initial population was divided into males and females.The two will use different update methods:

Table 5
Parameters of the solar PV plant.
Solar PV plant No.
Eq. ( 31) and Eq. ( 32) are the positional update formulas for male and female mayflies.

Movements of male mayflies
Considering that male mayflies are always performing a nuptial dance, hence they constantly small movements.pbest i is the individual with the best dominance level obtained so far.The velocity would be updated according to the dominance of x i and pbest i .The definition of dominance can be seen in Ref. [22].
The corresponding speed update formula is determined according to the dominance of the male mayfly in relation to the gbest, as expressed in the following Eqs.(33) and (34): where g denotes the progressively decreasing gravity coefficient.a 1 and a 2 are two different attraction constants.β denotes a visibility constant, d expresses the wedding dance coefficient, and Eq. ( 35) was used to calculate the Cartesian distance r p from x i to gbest.J. Zhu et al.

Movements of female mayflies
The choice of the velocity update formula for the i-th female mayfly y i is determined by the dominance relationship between y i and x i .Specifically, if y i is dominated by x i , the velocity is calculated as Eq. ( 36): where g is a gravity coefficient, r mf represents the Cartesian distance between x i and y i .Conversely, if y i dominates x i , y i would update the velocity by Eq. ( 37): where fl denotes the wandering coefficient.

Mating of mayflies
When female mayflies are attracted to males, they engage in mating behavior.According to random probability, male mayfly x i and female mayfly y i are selected to mate and generate two offspring mayflies (x off1 and x off2 ).The mating behavior can be formulated as Eqs.( 38) and ( 39): where G is a random value according to the Gauss distribution.

Performance of fast non-dominated sorting
According to the NSGA-II [22], a population is sorted into different non-dominated levels.Fig. 6 shows a hierarchical schematic diagram of non-dominated ranking.Specifically, each individual solution s comes with two additional elements: the total number of solutions dominated by this solution, denoted as n dom ; the set of solutions dominated by this solution, denoted as S p .Theoretically, we describe non-dominated sorting will perform uses O (MN 2 ) comparisons in the worst case.There are three main steps of the algorithm.
Step 1. Find all individuals s with n dom = 0 and save them in the set F1.
Step 2. For each individual s in the current set F1, find the set S p dominated by s, visit each individual q in the set S p and set n q = n q -1, if n q = 0, save the individual q in another set H.
Step 3. The first non-dominated front is F1 and set H as the initial set.
Repeat the above operations until the entire population being stratified.
J. Zhu et al.

Nonlinear decreasing inertia weight
In order to assist in achieving an adequate balance between exploration and exploitation, a nonlinear decreasing inertia weight w is introduced.In the early iterations of the NSMA, w tends to decrease slowly.The NSMA has a strong global search capability and can find better individuals faster.In the later iterations, w tends to decrease rapidly.At this time, the w value is smaller, which corresponds to a smaller search range of all the mayflies.The iterative update curve of the value of w is shown in Fig. 7.The specific calculation formula is displayed in the following formula (40): where w and w are the upper and lower bounds of w, iter and Max_iter represent the current number of iterations and the predetermined total number of iterations, respectively.

Superiority of feasible solutions (SF)
In [23], SF was mentioned for the first time.Michalewicz and Schoenauer have divided all common methods of constraint handling into five categories [24].SF is a combination of the 2nd and 3rd of these methods, which is a mixture of the penalty function method and the separation of infeasible solutions to obtain SF.In SF, the equality constraint is first processed, converting equality constraints to inequality constraints using Eq.(41).
whereθ identifies a tolerance parameter.H is the count sum of inequality constraints and G denote the number of total constraints.For any obtained feasible solution s i , s i must satisfy all inequality constraints C i (x,y), including those formed by the transformation of the equality constraints; otherwise, s i must be an infeasible solution.Next, an entity ε(x, y) is assigned to each infeasible solution s i to record the constraint violation degree of s i .The formula for ε(x, y) can be written as Eq. ( 42): where ϖ i (ϖ i = 1/C max,i (x, y)) represents a weight parameter.Finally, the decision to retain the s i is based on the objective function value of the s i and the constraint violation degree ε(x, y).
The flow chart of the NSMA-SF for MOOPF problem with stochastic renewable power is presented in Fig. 8.

Results and comparison
To verify the performance of NSMA-SF, the modified IEEE-30 and standard IEEE-57 bus test system were selected for testing, and ten different cases were conducted.MATPOWER, an open-source tools tool is used to simulate and optimize power systems.All code is Fig. 7.The iteration curve for w. done on MATLAB.In addition, the simulated tests have been run on a computer with Intel Core i7 CPU @3.4 GHz and 64 GB RAM.Each study case was run 20 times repeatedly.
The solution to the MOOPF problem is intend to obtain a best compromise solution (BCS) from the best Pareto front (PF) to minimize multiple different and conflicting objectives.The well-known hypervolume (HV) metric is used to access PF [24].In this experiment, the diversity of PF and the pros and cons of convergence are reflected by comparing the value of HV.The largest one is selected as the best PF, and a set of solutions is selected as the BCS among them.

Results and comparison in the modified IEEE-30 bus system
The system architecture diagram of 30-bus system is presented in the form of a single line diagram, which is shown in Fig. 2. As shown in the figure, it is worth mentioning that two wind-driven generators on buses 5 and 11, and a PV generator is fitted to bus 13.On the system, NSMA-SF is applied to resolve the MOOPF problem with stochastic renewable power.The modified 30-bus system parameter settings were confirmed before the experiment.
As aforesaid, the optimal results obtained by NSMA-SF on six separate cases are reported.First, different cases are analyzed and the values of the load bus voltage and control variables provided by the BCS are verified, and then the BCS calculated by the NSMA-SF algorithm is shown.The obtained BCS is compared with the results in other literature for the MOOPF problem.

Experimental parameter setting
A detailed description of the main structural components of the improved IEEE-30 bus test system is listed in Table 6.The allowable limit of each variable is presented in Table 7.These specified values above are consistent with those in other literature.The values for the best single target and BCS for all cases are also listed in Tables 8-9.It is obvious from Tables 8-9 and Fig. 9 that the constraint variables are all between the set upper and lower limits.The algorithm strictly adheres to the reactive power limit in all cases of the studies conducted in this paper.
During the test, the worst total VD values of the load bus are presented in Tables 8-9 Presumably, if the voltage limit can also be met with the maximum VD value, the load bus voltage will be safe for all other solutions with smaller VD values.At the highest VD, the distribution of voltage amplitudes on all buses (except for the generator bus) is shown in Fig. 10.It can be seen from the curves that NSMA-SF appropriately satisfies the constraints of the load bus voltage.

Comparative analysis of experiments
Tables 8 and 9 listed the final values of the control variables for the best solutions for a single objective and BCS for six study cases of the modified 30-bus system.From Tables 9 and it is obvious that all 11 control variables are limited within the standard range.
Table 10 reports the 20 run results of the HV values using NSMA-SF in cases 1-6 and compared with the HV values obtained from NSGA-II-SF [26] and ACNSDE-SF [26].The result comparison is based on the Max, Min, mean and Std among the 20 independent run times.In addition, nonparametric statistical analysis was performed using the Wilcoxon Signed Ranks test, where R + and R− refer to the sum of the rankings of NSMA-SF outperforming its competitors or its competitors outperforming NSMA-SF [27].According to Table 10, the NSMA-SF achieves better HV values in terms of Max, Min and mean compared with two competitive algorithms in all 6 cases.In terms of the Std, the values of NSMA-SF are smaller than those of NSGA-II-SF in cases 2-6, but there is no significant difference between them.Table 10 also lists the 20 runs of NSMA-SF and statistical results of the Wilcoxon rank sum test of two other competitive algorithms for HV values.From the Wilcoxon Signed Ranks test results, in test cases 1-6, the difference in HV values between the proposed NSMA-SF and the other two comparison algorithms is significant.To visually compare the HV value distributions of the three methods, the violin diagram is used to display the distributions.The violin diagram combines the advantages of box and density plots, it can display the distribution state and probability density of multiple sets of data.HV values distribution of three algorithms on cases 1-6 are in Fig. 11.In all 6 cases, the HV value assignment of NSMA-SF was significantly better than NSGA-II-SF and ACNSDE-SF.
It is worth mentioning that the optimal PF is determined by the maximum value of HV, the optimal PF obtained from three algorithms are shown in Figs.12-15.The choice of the BCS depends on the distribution of the non-dominated solutions along the PF.The quality of an optimal PF will directly affect the value of the BCS.As shown in Figs. 12 and 13, Case 1 and Case 2 are both 2-objectives optimization.According to the figure, the shape of PF is similar to a curve.It can be seen form Figs. 12 and 13 that the PF of NSMA-SF has a lower aggregation density.In Figs. 12 and 13, NSGA-II-SF and ACDSDE-SF tend to aggregate in local areas.NSGA-II-SF and ACDSEE-SF have poor coverage of true Pareto optimal fronts, so they cannot efficiently find as many global optimal solutions as possible.By comparison, it can be proved that NSMA-SF has better search performance in searching for the optimal PF.Therefore, it is fully proved that NSMA-SF is superior to the other two competitive algorithms in finding multiple Pareto optimal solutions and their corresponding optimal PF in terms of outcome measures in decision space and objective space.The NSMA-SF algorithm ultimately exhibits a very competitive search performance compared to the other two algorithms.
However, as the objective number increases, the difficulty of solving the problem increases exponentially, which results in a significant decrease in the performance of NSMA-SF.For the 3-objectives study cases, cases 3-5, whose true PF shape is a hypersurface, the distribution of PF can be shown in Fig. 14.In case 6, the number of optimization objectives is four, which is difficult to display with a graph.The parallel coordinate plot is used to display the obtained PF.The parallel coordinate plot is suitable for the analysis and comparison of multi-dimensional data (especially when the dimensions are greater than 3), and can directly observe the distribution of each variable.Fig. 15 is a parallel coordinate plot of the BCS of the best PF on case 6.The figure shows a reasonably good distribution of solutions for the four-objective optimization case study using NSMA-SF, and a wider distribution of PFs than those obtained with

Table 6
The structural components of the modified IEEE-30 bus system.

NSGA-II-SF and ACDSDE-SF.
In other words, the BCS obtained by NSMA-SF has a wider range of choices.In general, NSMA-SF can obtain a diverse and welldistributed PF for all cases.The distribution of PF in the decision space obtained by NSMA-SF is far superior to the results obtained by the other two algorithms.From the above discussion and analysis, the proposed NSMA-SF can provide good coverage of the real PF for the MOOPF problem with stochastic renewable power.
Due to the small amount of literature on integrating renewable energy generators in IEEE-30 bus systems, some results from the literature on solving classical MOOPF problems are selected for comparison.Table 9 presents the BCS for NSMA-SF and other optimization methods and compare the results.
In cases 3 and 4, compared with the BCS obtained by other algorithms, the emission in the BCS of MOOPF with renewable energy is the optimal (minimum).In the 3-objective optimization of case 5, the cost and loss values are the best, but with significantly high level of VD.In case 6, which is optimized with four objectives, the BCS obtained by the proposed NSMA-SF is better than other algorithms in optimizing the emission objective.In case 1, 3, 4 and 6, as long as emission is one of the optimization objectives, the emission of the BCS obtained by NSMA-SF is the best value.Overall, in all 6 cases, comparing the BCSs for solving the classical MOOPF problem, it can be found that the emissions in the BCSs obtained by the NSMA-SF are far smaller than their corresponding values.Explain that with the incorporation of renewable energy sources in the modified 30-bus test system can significantly reduce emissions.When comparing the optimal results of the simulation with those in other literature, it can be discovered that NSMA-SF can achieve better (smaller) emissions.It shows that NSMA-SF is an effective method for resolving the MOOPF problem with stochastic renewable power.J. Zhu et al.

Results and comparison in the standard IEEE-57 bus system
In order to verify the robustness of the NSMA-SF, a higher-dimensional standard IEEE-57 bus system was chosen to conduct simulation experiments.The details are given in Table 3.

Experimental parameter setting
Table 13 details the composition of the standard 57-bus system.In Table 14, the maximum and minimum values of all variables are presented.Fig. 16 presents the generator reactive power for all cases, showing that the algorithm complies with the generator reactive power constraints in all cases.Fig. 16 suggests that the reactive power of each generator is also strictly controlled within a preset range.The bus load for the four cases with the worst (maximum) cumulative VD of the load bus can be seen in Fig. 17.It can be found from the figure that the VD of each bus is also within the defined upper and lower boundaries.In addition, many load bus voltages were found to operate near the upper limit.This fact proves again that solving the constrained MOOPF problem requires an efficient constraint handling technique.

Experimental analysis in standard IEEE-57 bus system
Tables 15 and 16 list the final values of the control variables for the best solutions for a single objective and BCS for four case studies of the IEEE-57 bus system.Table 17 and Fig. 18 show the HV values and their distributions for cases 7-10, respectively.As the dimension increases, the HV value becomes smaller and smaller.On the other hand, the multi-objective solution increases    exponentially due to each additional number of objectives [42].With the increase in the number of objectives, the HV value decreases significantly, from 0.0648 to 0.0045.Fig. 19a and b are the optimal PF of cases 7 and 8 for 2-objective optimization, respectively.Fig. 19c is the optimal PF for case 8, where the number of optimization objectives is three.In this case, a PF with reasonably good diversity and distribution is achieved by the NSMA-SF algorithm.For the 4-objective optimization case 10 using NSMA-SF, a good distribution of Pareto solutions is also observed in the parallel coordinates plot in Fig. 19d.
In order to test and verify the ability of NSMA-SF in solving the MOOPF problem, we compare it with other multi-objective algorithms which are used to solve the same problem.Table 18 shows the comparison of the BCS obtained by NSMA-SF and the BCS obtained by other algorithms.It can be seen from Table 18 that no BCS can obtain a best value of every objective in multiple objectives than other algorithms.The BCS obtained by the algorithm does not represent the optimal solution of MOOPF, but it can provide manufacturers with a full reference scheme, so that they can make the final decision according to their own actual conditions.The BCS depends on the strategy and policy of the utility company based on its specific technical and economic forecasts.J. Zhu et al.

Conclusion and future work
In this paper, a constraint multi-objective algorithm, namely NSMA-SF is applied to resolve multi-objective OPF with renewable source.Generation cost, power loss, VD, and emission are selected as optimization objectives.The classic IEEE-30 bus system was modified by replacing the three thermal generators connected to it with two wind power plants and one solar photovoltaic power plant.To evaluate the performance of the proposed algorithm, the bus system is used to perform simulation experiments two different PDFs are used to calculate the actual power of uncertainty renewable energy sources.NSMA combined with SF has been successfully solved the MOOPF problem with stochastic wind and solar power.By checking the voltage of each load bus, this paper proves that the NSMA-SF algorithm can effectively.Experimental results show that the NSMA-SF provides a satisfactory BCS in different cases.Specifically, NSMA-SF achieves better PF diversity and Pareto Front (PF) convergence than NSGA-II-SF and ACDSDE-SF for the optimization cases of 2 and 3 objectives.This can also be verified from the distribution of PF.In the case of 4 objectives, NSMA-SF obtains a more uniform and wider distribution PF.From the comparison experiments, the replacement of conventional thermal generators with renewable energy generators has the potential to significantly reduce system pollution emissions and environmental pollution.It is verified that using NSMA-SF to solve the MOOPF problem is an effective approach on a higher dimensional IEEE-57 bus system.What should be pointed out is that the model presented in this paper considers the load of each bus to be deterministic.However, in reality, these bus loads are uncertain.In future work, we suggest that more constraint processing techniques can be used in MOOPF problem with wind and solar energy, such as self-adaptive penalty (SP), epsilon constraint (EC), stochastic ranking (SR) and tested on more complex bus systems.

Data availability
All data generated or analyzed during this study are included in this published article.

J
. Zhu et al.

Fig. 12 .
Fig. 12.The best PF obtained by three algorithms on Case 1.

Fig. 13 .
Fig. 13.The best PF obtained by three algorithms on Case 2.

Fig. 14 .
Fig. 14.The best PF obtained by three algorithms on Case 3.

Fig. 15 .
Fig. 15.The best PF obtained by three algorithms on Case 6.

Table 1
Cost coefficients of thermal generator in two test systems.

Table 2
Emission coefficients of thermal generator in two test systems.

Table 7
The maximum and minimum values for each variable in the modified IEEE-30 bus system.

Table 8
Simulation results obtained by NSMA-SF.
*Bold values represent the best values.

Table 9
Simulation results of BCSs in study cases 1-6.

Table 10
Statistical results of HV indicator value.
*Bold values represent the best results.

Table 11
Comparison of BCSs of cases 1-4 for IEEE-30 bus system.MOOPF cons b ders the valve point effect in the cost calculation.
a MOOPF with a renewable energy sources.b

Table 12
Comparison of BCSs of cases 5 and 6 for IEEE-30 bus system.MOOPF cons b ders the valve point effect in the cost calculation.
a MOOPF with a renewable energy sources.b

Table 13
The structural components of the standard IEEE-57 bus system.

Table 14
The maximum and minimum values for each variable in standard IEEE-57 bus system.

Table 16
Simulation results of BCSs in study cases 7-10.

Table 17
The value of the HV indicator obtained by NSMA-SF.
*Bold values represent the best results.bMOOPF cons b ders the valve point effect in the cost calculation.